Poisson Nernst-Planck Model of Ion Current Rectification 
through a Nanofluidic Diode 



O 

o 

(N 



D. ConstantiiE and Z. S. SiwjQ 

Department of Physics and Astronomy, University of California, Irvine, CA 92697 

(Dated: February 1, 2008) 

We have investigated ion current rectification properties of a recently prepared bipolar nanofluidic 
diode. This device is based on a single conically shaped nanopore in a polymer film whose pore walls 
contain a sharp boundary between positively and negatively charged regions. A semi-quantitative 
model that employs Poisson and Nernst-Plank equations predicts current-voltage curves as well as 
ionic concentrations and electric potential distributions in this system. We show that under certain 
conditions the rectification degree, defined as a ratio of currents recorded at the same voltage but 
opposite polarities, can reach values of over a 1000 at a voltage range < —2V,+2V >. The role 
of thickness and position of the transition zone on the ion current rectification is discussed as well. 
We also show that rectification degree scales with the applied voltage. 

PACS numbers: 81.05.Zx, 02.60.Cb, 81.07.De 
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I. INTRODUCTION 

Nanoporcs in polymer films and silicon materials have 
attracted a lot of scientific interest due to their applica- 
tion as biomimetic systems for models of biological chan- 
nels P, H, [1| , and as biosensors [1, [1, HI, 0] ■ Nanopores 
are also used as a basis for building ionic devices for 
controlling flow of ions and charged molecules in a so- 
lution 8]. In these systems electrostatic and hydropho- 
bic interactions between the ions and molecules that pass 
through the pore and the pore walls are used to amplify 
or stop the transport. Our group has recently prepared 
a nanofluidic diode [9] which rectifies ion current in a 
similar way as a bipolar semiconductor diode rectifies 
electron current. This diode is based on a single coni- 
cally shaped nanopore in a polymer film with openings 
of 2 -;- 6 nm and 1 fim, respectively. The surface charge 
of the pore is patterned so that two regions of the pore 
with positive and negative surface charges create a sharp 
barrier called the transition zone. This nanofiuidic diode 
is bipolar in character since both positively and nega- 
tively charged ions contribute to the measured current. 
Another type of nanofiuidic diode was prepared in the 
group of A. Majumdar [10] with a sharp barrier between 
a positively charged zone and a neutral part of the pore. 
Presence of only one type of surface charges causes the 
latter device to be unipolar. 

In this article we provide a semi-quantitative descrip- 
tion of ion current rectification of the bipolar diode 
based on single conical nanopores. We use contin- 
uum mean-field theory, and more precisely Poisson and 
Nernst-Planck (PNP) equations, to define our model 
[ill [T^ . |T3| that we apply to describe current-voltage 
characteristics as well as profiles of ionic concentrations 
and electric potential in the system. Our model is 
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based on a similar approach as described by Cervera 
et al. in [l^, [l^. We would like to mention that 
the model of nanofiuidic diode based on a rectangular 
channel with limiting length of tens of nanometers pre- 
sented in [l3| combines PNP with Navier- Stokes equa- 
tions. Nanopores considered here are significantly nar- 
rower with the smaller diameter around 5 nm. There- 
fore, we consider here only PNP equations, because in 
nanopores of diameters less than 10 nm, electrophoresis 
is known to be the most dominant process [13, [l!] ■ 

It is important to mention that continuum description 
of ion transport by the PNP equations has successfully 
rendered main properties of many biological channels, 
e.g. see [l^ for the calcium channels treatment. When 
compared to more precise but computational expensive 
techniques, e.g. Brownian dynamics [20| . the continuum 
approach has a clear advantage in terms of computing 
cycles. As discussed in [2^ this description produces 
results which correctly render the physical and chemi- 
cal phenomena occurring inside a cylindrical nanopore 
whose diameter is bigger than 2 Debye lengths of the 
considered ionic species. Hence it is plausible to obtain 
valid results for the conical nanopores as long as their 
smallest diameter is bigger than 2 Debye lengths in the 
examined system. 

In this study, we thoroughly discuss how ion current 
rectification in a nanofiuidic diode depends on the po- 
sition of the transition zone between the two regions of 
pore walls with opposite surface charges. Infiuence of 
the width of the transition zone on rectifying properties 
of the device is studied in detail as well. We identify the 
range of values for the position of the transition zone and 
its width that lead to a significant enhancement of the 
current rectification degree. This aspect is of crucial im- 
portance for further improvement of these new devices. 
The paper is organized as follows. In section|TT|we explain 
the experimental system and present the theoretical basis 
of our model. In section IIIII we present the results of the 
analysis and provide directions to future improvements 
of the design of nanofluidic diodes. 
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II. METHODS 

In this section we present the experimental results and 
the theoretical background which provide the basis of 
our simulations. In the first subsection we briefly explain 
the experimental facts which inspired the current project. 
The following subsections are devoted to the theoretical 
background. 



A. Experimental 

Single nanopores in polymer (polyethylene terephtha- 
late, in short PET) film were prepared by the track- 
etching technique as described in [2l|, Briefly, the 
method consists of irradiation of polymer films with sin- 
gle swift heavy ions [1^ (Gesellschaft fuer Schwerionen- 
forschung, Darmstadt, Germany), and subsequent asym- 
metric etching of the irradiated foils at 9M NaOH 21]. 
This procedure leads to preparation of membranes that 
contain only one pore with the small opening, called tip, 
as narrow as several nanometers in diameter, and the 
big opening, called base, with micrometer size diame- 
ter. As the result of the chemical etching, there are car- 
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FIG. 1: Pore geometry with schematic representation of sur- 
face charge distribution creating a bipolar nanofluidic diode. 

boxyl groups formed on the surface of pore walls with 
density of about one per nrri^. which determines the 
surface charge to — 0.16 C/m^ [2J]. In order to mod- 
ify the surface of PET nanopores we used the method of 
coupling carboxyl and amine groups with l-Ethyl-3-[3- 
dimethylaminopropyljcarbodiimide hydrochloride (EDC) 
[25(1 . The reaction between the carboxyl groups and di- 
amines changes the surface charge from negative to pos- 
itive. Targeted modification of the tip was obtained 
by introducing the reagent mixture only on the side of 
the membrane with the small opening Q. The other 
side of the membrane was in contact with a buffer so- 
lution. A conical shape of the pore makes the distri- 
bution of the diamines and EDC in the pore extremely 
non-homogeneous. There is a high concentration of the 



reagents at the tip of the pore that assures the chemical 
modification reaction of the tip to occur. The concentra- 
tion of diamines and EDC at the remaining part of the 
pore is very low, therefore this part remains negatively 
charged. As the result, we obtain a surface charge dis- 
tribution as schematically shown in Fig. [T] An example 
of the current-voltage characteristic through such a de- 
vice recorded at 0.1 M KCl and pH5.5 is presented in 
Fig. [2l The rectification degree of the system is defined 
as the ratio of currents recorded for positive and negative 
voltages, respectively 



f{V) = 



\I(V)\ 

\ii-v)\ 



(1) 



The rectification degree of this system (Fig. ^ equals 
217 at 5V [si]. We would like to mention that even de- 
vices with uniform surface charge distribution rectify the 
current but their rectification degree is one order of mag- 
nitude smaller than the one obtained for the systems de- 
scribed above 0, [2l|. It is also important to mention 




Voltage [V] 

FIG. 2: Experimental data of ion current recorded at 0.1 M 
KCl, pH 5.5 through a single conical nanopore with surface 
charge distribution as shown in Fig[T][||. The pore diameters 
were 5nm and 1000 nm, respectively 

that the single pore diode is a nano-realization of bipolar 
membranes consisting of cationic and anionic membranes 



B. Modeling 

Before we start to introduce the theoretical base of the 
model we would like to specify the notation. We have 
assigned the subscripts "T" and "i?" to designate the 
various physical quantities computed at the tip and at 
the base sides of the pore, respectively. In this sense 
is the location of the pore tip and is the location of 
the pore base. The thickness of the membrane is denoted 
with L. The cone opening angle is 9m- The membrane 
is immersed in an electrolyte containing potassium and 
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chloride ions, and in addition we apply an external elec- 
tric field with the help of two electrodes positioned far 
away from the membrane on each side of the pore. Due to 
Gn, resistance of single nanopores, the voltage drop will 
occur only inside the pore, not in the electrolyte bulk. By 
convention, we will keep the electrode on the pore tip side 
grounded such that the sign of the electric potential dif- 
ference will be controlled by the electrode on the base side 
of the pore. The existence of the external field causes the 
ions to migrate towards the electrodes, which we record 
as ion current. We would like to mention that for small 
pores [2| and more complex electrolytes in which chemical 
reactions can occur, one observes large fluctuations of ion 
current in time [HjllO]- However our model in its present 
form docs not include chemical reactions between com- 
ponents of the electrolyte, and therefore the system will 
reach a stationary regime, i.e., the ionic fluxes are con- 
stant in time. In nanofluidic diodes, constant currents 
at constant voltages are indeed observed experimentally. 
Our goal is to compute the steady state current for a 
fixed set of physical and chemical parameters. Hence we 
express the stationarity of the ionic fluxes using Nernst- 
Plank condition. To be more specific we assume that for 
each ionic species, indexed by "i" , the molar flux deflned 
as 

J, (f) = - A [ Vc, (f) + (f) V0(f) ] , (2) 
is conserved, i.e., it obeys the Nernst-Plank condition 



VJ,(f) = 0, 



(3) 



where Dt is the diffusion coefficient of an ion i, a is the 
molar concentration and Zi is the charge number of the 
ion i. We have denoted with cj) the electric potential in 
RT/F units, which satisfies the Poisson equation 



A(/)(f) = 



(4) 



where e is the electric permittivity of the medium and 
the summation is carried over all the ions present in the 
solution. The parameter F is the Faraday constant, R is 
the gas constant and T = 293 K is the absolute tempera- 
ture. The PNP equations ^ and ^ provide a complete 
set of equations which together with appropriate bound- 
ary conditions allow us to determine the concentration 
of each ionic species and the electric potential inside the 
pore. Consequently the total ion current density along 
an arbitrary direction h can be computed with the help 
of © 



J 



(5) 



where the summation is performed again over all consid- 
ered ions. 

As the problem has azimuthal symmetry with re- 
spect to the pore axis, it gets reduced to a set of two- 
dimensional second order partial differential equations. 



To simplify the problem even more we perform an ana- 
lytical integration with respect to the polar angle and we 
obtain a set of second order ordinary differential equa- 
tions which is finally solved numerically. We would like 
to mention that after integration one has to neglect one 
term in order to obtain a well defined set of ordinary dif- 
ferential equations ^16^. Our comparison tests between 
the full two dimensional problem and the approximated 
one dimensional problem show insignificant differences in 
both the quantitative and the qualitative results. Hence 
the solution of the problem will provide averaged quan- 
tities in the sense that the polar angle dependence is lost 
during the integration. However our main interest is to 
compute the ion current, which involves averaging along 
the polar coordinate, therefore we do not loose any in- 
formation. To be more specific, the differential equations 
which result from the above described procedure have the 
following structure 



d r 
dr 



dci{r) d(j){r) 
+ ZiCi{r) 



and 



1 d 



dr 



d4){r) 



dr 







(6) 



dr V dr 



eRT 



[Y,z^c,{r)+X{r)], (7) 



where the quantities Ci{r) and (/)(r) in ^ and ([7]) come 
from integration, i.e.. 



1 



1 - cos 9m Jo 
1 



1 - cos Om Jo 
and the function X{r) is given by 
sin (?Afc(r) 



dO sine (j){r,t 
dO sin Ci (r. 



X{r) = 



Fr{l-cos0M) ' 



(8a) 
(8b) 

(9) 



where a{r) represents the surface charge density. For 
simplicity we have kept the same notation for the aver- 
aged quantities Ci(r) and 0(r) and we will use the func- 
tional dependence on the polar angle to distinguish be- 
tween them and the non-averaged quantities Ci{r, 6) and 



C. Boundary Conditions 

To completely define the problem we have to spec- 
ify the boundary conditions. These conditions are de- 
rived from the so called Donnan equilibrium conditions 
which express the thermodynamic equilibrium between 
the membrane with excess surface charge and the ionic 
concentrations in the bulk. One can derive these condi- 
tions as in [sH and obtain the following relations for the 
concentrations at the pore border 



Ci{rj) = coi exp[-2;i0£)(rj)] 



(10) 
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where Coi represents the bulk concentration of ion i and 
j — T,B indexes the tip and the base sides of the pore, 
respectively. (j>D is the Donnan potential defined as the 
potential difference created across an ion exchange mem- 
brane 



(11) 



'^z,Ci{rj) + X{rj) = 0, 



(12) 



provides an equation for exp (rj )]. Hence, the bound- 
ary conditions can be derived from (fTU]) and as 

c,{r,)^^[-z.Xirj} + ^/x{^^^y7^,), (13) 



and 



Zi Co 



(14) 



The PNP equations ([6]) and ([7]) together with the bound- 
ary conditions ([T5)) and HH) form the base of our model. 
We deal with two ionic species in our system, i.e., 
and Cl~ , therefore we have to solve a set of three second 
order ordinary differential equations. 



D. Surface Charge Distribution 

The charge distribution is formed by two distinct re- 
gions with positive and negative surface charges respec- 
tively, and a transition zone between them. These regions 
are characterized by the asymptotic values of the surface 
charge density, i.e., the value of a far away from the tran- 
sition zone. We employ some simplifying assumptions in 
order to characterize the surface charge distribution in 
the transition zone. Firstly, we consider that the magni- 
tudes of the asymptotic values of the positive and neg- 
ative surface charge densities are the same and equal to 
eo/nrn^, where cq is the elementary charge. Additionally, 
we assume that the transition between the positive and 
negative charges on pore walls is made in a symmetric, 
continuous, and monotonic way. Hence we will consider 
the surface charge density in the pore to be described by 
the logistic type function (see Fig. [3]): 



l + cxp[-2fc(r-ro)] 



(15) 



where ctq — —eo/nrn^ is the minimum value of the sur- 
face charge density, rg is the position along the pore axis 
where cr(r) equals zero, and k measures the slope of a in 
the transition zone. Another equivalent way to charac- 
terize the charge distribution in the transition zone is to 
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where (j>j is the electric potential at the corresponding o -0.05 
electrode. The electro-neutrality condition, which has 
the form 



U -0.10 



-0.15 



w=20nra 



=1.0e-06 



_L 



9.00e-07 9.50e-07 l.OOe-06 1.05e-06 
Position Along Nanopore Axis [m] 



l.lOe-06 



FIG. 3: Surface charge distribution that has been used in 
calculations. The width of the transition zone w and position 
of the zero charge point ro are indicated in the figure. 



specify the transition zone width w (see Fig. [3]) which is 
defined as the part of the pore where the magnitude of 
the charge is bellow |tTo|/e, i.e.. 



— In — 

2k mi 



(16) 



where mi^2 = tanh^^(l/2) are constants. In this problem 
there are many parameters which can be modified. How- 
ever our main interest was to see how the properties of 
the nanofluidic diode modify when the surface charge dis- 
tribution parameters change. Therefore, we have chosen 
to keep all the parameters constant, the only exception 
being rg and w, the parameters, which we believe can be 
controlled experimentally. 



III. RESULTS 

As we have discussed in the previous section we have 
to solve a boundary value problem (EVP) defined by a 
set of three second order ordinary differential equations. 
The solution will provide the potassium and chloride con- 
centrations and the electric potential as functions of the 
coordinate along the pore axis. We have used the EVP 
solver described in . 32i] to compute the solution with an 
error less than 10~^. The solution of the computation is 
used to evaluate the ion current which is given by 



27r(l 



(17) 



where jr is the radial component of the ion current den- 
sity given by (ini) . As we expect the current is not depen- 
dent on the radial coordinate and it has a constant value 
along the pore axis. 
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FIG. 4: Profiles of ionic concentrations and electric potential in a nanofluidic diode based on a single conical nanopore for 
To = 1 ^im, w — 100 nm, and various voltages as indicated, (a) profiles of Cl~ concentration in a nanopore for forward bias 
(positive voltages); (b) concentration profiles for forward bias; (c) profiles of Cl~ and K'^ concentrations for reverse bias 
at —IV; (d) electric potential profiles in the nanopore for forward and reverse bias. 



A. Current-voltage curves and ionic concentration 
profiles 

The system that we have considered is a single conical 
nanopore with pore diameters 5 nm and 1000 nm, respec- 
tively, at symmetric electrolyte conditions of 100 mM 
KCl on both sides of the membrane. The thickness of 
the membrane is 12 ^m. Fig. IH^a-c) show the numerical 
solutions of the corresponding distributions of potassium 
and chloride ions in forward and reverse bias, respec- 
tively. Fig. md) presents the distribution of the electric 
potential in the nanopore for various applied voltages. 
The increase of ionic concentrations in the forward bias 
is remarkable. Concentrations of both ions, potassium 
and chloride, have increased by an order of magnitude, 
compared to the bulk concentration of 100 mM. This re- 
sult also confirms bipolar character of our device, since 
and Cl~ contribute almost equally to the measured 
ion current. Ionic concentrations at the reverse bias are 



dramatically different. We see significantly lower concen- 
trations of both ions in a region close to the transition 
zone which indicates the formation of the depletion zone. 
The inset in Fig. [4fc) shows a magnified view of the de- 
pletion zone where the concentration of both ions drops 
practically to zero. This region is located slightly off the 
center of the transition zone because of the conical geom- 
etry of the nanopore. Fig. HJ^d) shows profiles of electric 
potential in the pore for various voltages applied in the 
forward and reverse directions. Similar to the results for 
diodes based on rectangular channels and bipolar mem- 
branes [l3, [H, [13, HI], in our bipolar diode based on 
single conical nanopores, the whole voltage drop occurs 
in the depletion zone. This sudden drop of the electric 
potential causes the appearance of a huge electric field 
located in the depletion zone, which will effectively block 
the ionic currents in the reverse bias. 

Using the numerical solution and Eq. p7p we com- 
pute the ion current. Fig. [5] shows a typical current- 
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FIG. 5: Current-voltage curve for a nanofluidic diode calcu- 
lated with Eq. (fT7)l for ro = 100 nm and w = 100 nm. The 
surface charge was assumed as shown in Fig. O The pore 
openings were 5 nm and 1000 nm, respectively. 



voltage curve calculated for this nanoporc with param- 
eters To — 100 nm and w — 100 nm. Similar to the ex- 
perimental data shown in Fig.[2l the model predicts very 
high rectification degrees with very small ion currents for 
reverse bias. 



Position of depletion zone, dynamics of 
nanofluidic diode formation 



The process of preparation of nanofluidic diode is based 
on a chemical reaction of the reagent with the carboxyls 
groups on the pore walls. This causes change of the po- 
sition rg and increase of the rectification degree in time 
when the chemical reaction progresses [9]. We decided 



l.Oe-07 



< 

S 5.0e-08 



O.Oe+00 











— r|j=10nm 






— r|j=50nm 

— r|j=75nm 
i|=100nm 

ijj=160nm 






rg=1000nm 














-4 





Voltage [V] 



FIG. 6: Current-voltage curves for a nanofluidic diode (Fig. 
m for various values of the ro parameter when w = 100 nm. 
Note the change of the direction of ion current rectification. 

to investigate in more details how sensitive the rectifica- 



tion degree is on the position of the transition zone in 
the pore. We performed therefore numerically the exper- 
iment of chemical modification of the pore walls. How- 
ever, the reader must be aware that our surface charge 
distribution is not necessarily the distribution which cor- 
responds to the real case. In the real case the shape of 
cr(r) changes in time i.e., it changes with tq, whereas 
in our case, the time evolution is simply described by a 
translation of a{r) along the radial coordinate. We think 
however that the main features of the real charge distri- 
bution are present in our model system and this gives 
us the possibility to understand and get insight into the 
nanofiuidic diode behavior when the surface charge gets 
modified. A given value of the transition zone width w 
was assumed and the position rg was moved between the 
pore tip and its base. In other words, we start from the 
situation where the whole surface of the pore is negatively 
charged and finish in a state where the whole surface is 
positively charged. In the initial and final states we have 
a homogeneously charged conical geometry that leads to 
ion current rectification via rocking ratchet mechanism 
described earlier in [s^ and [s^l . The rectification degree 
was defined in Eq. ([T]) as the ratio of currents recorded 
for positive voltages divided by currents recorded for neg- 
ative voltages. These homogeneously charged conical 
nanopores are characterized by rectification degrees less 
than 10. As we move the position tq inside the pore the 
surface charge in the tip changes acquiring more positive 
charges. As shown in Fig. [SI it is remarkable that a small 
amount of positive charge in the tip of the pore reverses 
the direction of rectification and changes the devices from 
monopolar to bipolar. The rectification degree versus rg 
is shown in Fig. [7] for several values of the transition 
zone width. The rectification degree increases steeply as 
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FIG. 7: Rectification degree of a nanofluidic diode as a func- 
tion of "ro" parameter for various thicknesses of the transition 
zone 'W (Fig. El, at 2V. 

we move rg inside the pore along its axis, reaches a max- 
imum, and decreases to a value lower than 10, when rg 
is located at the base of the pore. A very strong depen- 
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dence of the rectification degree on tq, when tq is located 
at the narrowest part of the pore, suggests that transport 
properties of the system are determined by the physical 
and chemical properties of the tip of the nanoporc, in this 
case the surface charge distribution. It is also important 
to mention that we have predicted higher currents for 
devices where vq is located further from the tip (e.g. at 
rg = lfj,m). However, these devices show lower rectifi- 
cation degrees due to larger leakage currents (off state 
currents). Depending on the specifics of the various ap- 
plications nanofluidic diodes can therefore provide higher 
(lower) currents and lower (higher) rectification degrees. 
We would like to emphasize that all these various prop- 
erties can be obtained by just changing the position of 
the zero surface charge point. 

As a comparison, we have also included rectification 
data for a monopolar diode designed such that the tip 
of the pore was positively charged (0.16 C/m^) with the 
rest of the pore neutral (OC/m^). The surface charge 
distribution is similar to the one shown in Fig. [31 the 
only difference is that tq represents in this case the posi- 
tion where the surface charge is half its maximum value. 
The parameter w was set to 100 nm. One can see that 
overall the rectification degrees for the monopolar diode 
are much smaller than in case of the bipolar diode when 
both devices are studied in the same KCl concentration. 

Our results show that a huge increase of ion current 
rectification degrees is possible when the surface of a pore 
is patterned such that a part of a pore with negative sur- 
face charge is brought into contact with a part of the pore 
with positive surface charge. Our analysis also points to 
the importance of future experimental efforts aimed at 
a better control over the width of the transition zone as 
well as the position of the center of a depletion zone. 

C. Scaling of rectification degree with voltage 

We have also analyzed the behavior of the rectification 
degree with voltage. In Fig. [8] we show the rectification 
degree of a nanofiuidic diode f{V) as a function of voltage 
ranging between —5 V and 4-5 V for the case when the 
transition zone is located 1 fim away from the tip, and 
for various values of w. We find that f{V) has a power 
law dependence for voltages larger than « 0.3 V, 

f{V)^V\ (18) 

where the scaling exponent 7 w 1.81. For guidance, the 
power law with exponent equal to 2 is shown in Fig. [S] 
The power law scaling of the rectification degree with 
voltage might suggest the existence of universal scaling. 
To our knowledge, this is the first time when this scaling 
has been observed in the variation of the rectification 
degree with the applied voltage for a nanofiuidic diode 
based on a single conical nanopore. The applicability 
of the universal scaling and universal scaling exponents 
to nanofluidic diodes will be investigated in our future 
studies. 
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FIG. 8: Log-log plot of the rectification degree of a nanofluidic 
diode (Fig. [TJ as a function of applied voltage at ro — Ifim 
and for various values of the "w" parameter. The red dashed 
line is a power law with exponent equal to 2. 



D. Dependence of depletion zone width on voltage 

Another interesting question to consider is the depen- 
dence of the depletion zone width Wd on the reverse bias 
voltages. The depletion zone is defined as the region 
where mobile charges are absent or at least their number 
is small compared to adjacent regions Fig. [H This ion 
free region appears only in the reverse bias regime. The 
lack of ions in this zone produces a drop of the electric 
potential, or in other words the electric field in this re- 
gion is strong enough to block the migration of mobile 
charges through the pore. To help visualize the above 
statements we have superimposed in Fig. |9] the concen- 
tration profiles, the electric potential, and its first deriva- 
tive which is proportional to the electric field. Since the 
voltage drops over a very small distance the electric fields 
created are of the order of magnitude lO^V/m . As we 
have explained above, the electric field has a strong peak 
in the depletion zone and we have considered the width 
of this peak to provide the definition for the width of the 
depletion zone. In order to plot all these profiles in one 
figure, scaling of their values was necessary. The electric 
potential was scaled up 50 times and its first derivative 
was scaled down 5000 times such that all the graphs will 
provide sufficient details. The plots in Fig. [5] correspond 
to w = 100 nm, tq = 250 nm, and for a reverse bias volt- 
age of —2 V. 

We would like to emphasize the difference between the 
width of the transition zone w and the width of the deple- 
tion zone. Larger w values lead to larger values of Wd, as 
shown in Fig. 1101 For w equal to 1 nm and 10 nm, the re- 
sulting Wd for 5 V exceeds the value of the transition zone 
width that assured formation of much smaller leakage 
currents, and consequently higher rectification degrees. 
For transition zone of width 100 nm and 200 nm, the de- 
pletion zone is narrower than the values of w. Larger Wd 
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FIG. 9: Profiles of ionic concentrations, electric potential and 
its first derivative at the depletion zone created at the reverse 
bias of —2 V. Values have been scaled (see text) to make the 
graphic visualization possible. 
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FIG. 10: Dependence of depletion zone width Wd on voltage 
for various values of the width of the transition zone w. 

therefore does not necessarily imply higher rectification 
degrees, as explained in the previous section. Smaller 
values of w and Wd assure larger electric fields for the re- 
verse bias, which also helps formation of a depletion zone 



of a larger resistance. 



IV. CONCLUSIONS 

The modeling of a nanofluidic diode presented in this 
manuscript has been motivated and inspired by recent 
experimental realization of this device [9|. In this arti- 
cle we have shown that rectification degrees of 10^ can 
be achieved with this device upon more precise surface 
charge control. We have shown that the width of the 
transition zone and the position of the transition point 
are parameters which influence very strongly the recti- 
fication properties of nanofluidic diodes. These param- 
eters in turn are controlled by chemical modiflcation of 
pore walls. Our analysis also shows that preparation of 
a nanofluidic diode based on a conical geometry offers 
much freedom in tuning rectification properties of the 
device. Future experimental efforts will be focused on 
developing surface chemistries so that various rectifica- 
tion degrees can be obtained. Our modeling also provide 
directions for the future design of the nanofluidic diode 
in order to prepare versatile devices which meet the re- 
quirements of various applications. 

It would be very interesting to see how the rectifica- 
tion behaves when other parameters of the problem are 
changed, for example shape of the nanopore, and magni- 
tude of surface charge. Another interesting aspect of the 
bipolar diode is the power law dependence of rectification 
degrees on voltage, which might suggest existence of uni- 
versal scaling and critical exponents in the system. Our 
future studies will be directed towards identification of 
these critical exponents and relating them to formation 
of the depletion zone. 
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